Nonadiabatic effects in the dynamics of atoms confined in a cylindric 
time-orbiting-potential magnetic trap 
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In a time-orbiting-potential magnetic trap the neutral atoms are confined by means of an inho- 
mogeneous magnetic field superimposed to an uniform rotating one. We perform an analytic study 
of the atomic motion by taking into account the nonadiabatic effects arising from the spin dynam- 
ics about the local magnetic field. Geometric-like magnetic-fields determined by the Berry's phase 
appear within the quantum description. The application of a variational procedure on the original 
quantum equation leads to a set of dynamical evolution equations for the quantum average value of 
the position operator and of the spin variables. Within this approximation we derive the quantum- 
mechanical ground state configuration matching the classical adiabatic solution and perform some 
numerical simulations. 

PACS numbers: 32.80.Pj, 03.65.Bz, 45.50.-j 



I. INTRODUCTION 



The difficulty of analysing a complex physical system 
is greatly reduced when one is able to identify a cer- 
tain number of different time scales present in the sys- 
tem dynamical evolution. Thus, a series of approxima- 
tions, generically termed adiabatic approximations, can 
successfully be carried out. The very simple basic idea 
is that of dealing first with the motion of the fast vari- 
ables, keeping the slow ones fixed but arbitrary, and 
then to complete the analysis of the entire system by 
allowing a variation of the previously fixed coordinates. 
The quantum adiabatic theorem and the molecular Born- 
Oppenheimer approximation are well-known examples of 
this approach, with its origins in the early days of quan- 
tum mechanics. The quantum-adiabatic theorem dic- 
tates that a system prepared in an eigenstate of its Hamil- 
tonian will remain in the corresponding eigenstate as the 
Hamiltonian is varied slowly enough. If the Hamiltonian 
returns to its original form, the system assumes the orig- 
inal eigenstate multiplied by an appropriate dynamical 
phase factor related to the instantaneous eigenvalue of 
the Hamiltonian. Berry made the interesting observa- 
tion that in addition to the dynamical phase factor pro- 
duced by the eigenvalue time evolution, the wavefunction 
acquires an additional phase contribution 1]. This ad- 
ditional contribution, the geometric phase, depends only 
on the path travelled by the system in the space of ex- 
ternal parameters. 

A canonical example for a system where this behav- 
ior occurs is that of a neutral particle carrying a mag- 
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netic moment and moving in an inhomogeneous magnetic 
field. Here, the fast variable is the transverse magnetic 
moment, and the slow variable is the atomic position 
and momentum 0. If the magnetic field varies slowly 
enough in space, the effective Hamiltonian governing the 
dynamics of the slow external variables contains an in- 
duced gauge potential, the so-called geometric potential. 
In the classical limit the gauge geometric fields acting on 
the neutral particle with a magnetic moment have been 
studied by Aharonov and Stern they found that the 
atom experiences geometric Lorentz-type and electric- 
type forces 0. The magnitudes of these forces do not 
depend on the amplitude of the magnetic field, but only 
on its local orientation. 

In order to treat the non-adiabatic corrections im- 
proved Born-Oppenhcimer methods were introduced for 
the case of arbitrary spin values Later the non- 

adiabatic terms modifying the atomic motion have been 
studied by several authors in the context of magnetic 
structures guiding or confining very cold atoms |(J 0, 
0, El- The non-adiabatic corrections produce spin- flip 
transitions leading to atomic loss from the magnetic con- 
figurations, and also they modify the atomic motion. 
High order post-adiabatic corrections, leading to geo- 
metric electromagnetism potentials, have been investi- 
gated for the elegant configuration of an atom orbit- 
ing around a straight current-carrying wire [Io|. For 
the verification of Berry's phase and its consequences 
a natural question is whether one can observe the di- 
rect modification of the atomic motion in the classical 
limit for the induced gauge potentials. Measurements on 
the motion of a rubidium Bose-Einstein condensate in 
a time-orbiting-potential (TOP) magnetic trap represent 
a quite stron g in dication for the existence of these geo- 
metric forces [TJ D3 ■ Those observations were analyzed 
through a classical description for the condensate center- 
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of-mass motion and for the atomic magnetic moment. 

In the present paper we perform an analytic study 
of the atomic quantum dynamics within TOP magnetic 
traps. We take into account the non-adiabatic effects 
arising from the dynamics of the spin orientation around 
the local magnetic field. Within a pure quantum descrip- 
tion, the geometric magnetic fields appear as a conse- 
quence of the presence of inhomogeneous magnetic fields. 
In this context spinors quantities can be introduced to 
describe the atomic spin states, as done by Ho and 
Shenoy |(| for the Berry phase in atomic condensates in 
magnetic traps. We derive an effective atomic dynamics, 
by means of a time-dependent variation principle making 
the quantum description analytically treatable. Thus, 
the atomic motion results by the coupling of a quan- 
tum harmonic motion, governing the atomic scalar wave 
function, and an effective nonlinear spin dynamics. The 
harmonic Hamiltonian depends on time-varying param- 
eters that, in turns, are linked to the spin state. Also 
the spin Hamiltonian parameters are time-dependent, as 
they result from the atomic wave-functions expectation 
values of the above geometric operators. Within this non- 
adiabatic approximation the ground state configuration 
matches well the adiabatic solution. We also performed 
numerical simulations with these new equations. For the 
parameters suggested by the standard experimental set- 
up we found that the adiabatic approximation is well 
suited. Nevertheless, by reducing the intensity of the 
bias field, non-adiabatic effects show up, because under 
these conditions the influence of the geometric fields is 
more relevant. 

Section [H] summarizes the classical analysis based on 
the adiabatic approximation and leading to the atomic 
micromotion. Section II I II reports a quantum analysis 
of the atomic motion within the adiabatic approxima- 
tion by taking into account the lowest frequency terms 
of the time-dependent potential. Within this approxima- 
tion we recover the quantum counterpart of the classical 
micromotion. Section IIVI studies the quantum dynam- 
ics of atoms into TOP traps. In Section by means of 
a timc-dependent-variation-principle, we derive an effec- 
tive dynamics for the atomic motion. Section IVII reports 
numerical simulations for the dynamical regime of the 
atomic motion. 



II. TOP TRAP 

Bose-Einstein condensation in dilute atomic gas is cre- 
ated by trapping cold atoms in a magnetic trap of which 
the Ioffe-Pritchard (IP) and the time-orbiting-potential 
are the most common ones. In a TOP trap the magnetic 
field, schematically represented in the inset of Fig. ^ 
is composed of a quadrupole (inhomogeneous) field and 
a rotating (time-dependent) bias field, Bo. The TOP 
trap, introduced in |13| for the very first experiments on 
Bose-Einstein condensation |14|. is employed by a num- 
ber of research groups producing Bose-Einstein conden- 



sates [II EE El El 113 

The single particle Hamiltonian for the atoms inside a 
magnetic field configuration which characterizes the trap 
geometry is given by 

H(t) = ^+mgz-^s-B(K,t) (1) 
2m s 

where s are the spin operators of the s = hj represen- 
tation, the last term takes into account the magnetic in- 
teraction energy of an atom with magnetic moment fis/s 
and projection /i = — along the magnetic field B. We 
also adopt the representation x — xx + yy + zz for the 
position vector. For a TOP trap the magnetic field is 
the superposition of a static quadrupole field and one 
rotating at the radiofrcqucncy (RF) lot 

B(x,t)=b(x)+B t (i). (2) 
Its components are 

b(x) = b x xx + b y yy + b z zz . (3) 

and 

Bi(t) = B cos(uiTt)x + B sm(uiTt)y ■ (4) 

The magnetic field parameters define the specific type of 
TOP trap we are analyzing 21]. By supposing the RF 
field rotating in the horizontal x,y plane we define the 
TOP geometry of the traps operating at Boulder ^| and 
at Pisa [22j. In this work we analyze the dynamics of a 
cylindric TOP trap with 

b X =by = -b Z /2 = b. (5) 

For the usual TOP trap three different time scales ex- 
ists: the fastest motion, given by the Larmor frequency 
LOjj and related to the spin precession around the lo- 
cal magnetic field; the magnetic bias field rotating fre- 
quency ljt giving rise to time-dependent forces at the 
frequences u)x, 2lot 1 ■ ■ ■ and a slower motion associated to 
the atom spatial motion given by the trap harmonic fre- 
quency LOh- These time scales are in order of magnitude, 
i.e. ujh << ljt << ljl thus making possible the adiabatic 
approximation. In the first place the fast spin precession 
around the local magnetic field allows to consider the 
atoms spin locked to the local magnetic field throughout 
the whole spatial motion. This leads to an adiabatic time 
dependent potential U = — ■ B(x,i) = |//| • |B(x,t)|. 
Secondly by averaging in time (over a period 2tt/lot) 
this potential and by keeping only the slowest component 
gives rise to a harmonic potential spatially confining the 
condensate. The next order of approximation examines 
the fast variables (fast with respect to the harmonic dy- 
namics at frequency u>h) related to the time-dependent 
potential (at frequency lot)- An exhaustive computation 
of this approximation is found in 23] . 

Gov et al [2(| have used the standard classical equa- 
tions of the atomic motion of a magnetic moment within a 
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Anti-Helmholtz coil 



micromotion orbit 




Anti-Helmholtz coil 



FIG. 1: Schematic representation of a cylindric TOP trap 
with the anti-Helmholtz coils producing a quadrupole field 
with vertical symmetry, the -Bo bias field rotating at angular 
frequency Wt in the horizontal plane, and the atomic cloud 
(black sphere) following the micromotion orbit. 



TOP in which all the different time scales are present. In 
this context a steady periodic orbit can be found exactly 
without resorting to any approximation. This periodic 
solution corresponds to what it is known as the atomic 
micromotion |2J] (see Fig. . As the atoms trace out 
the atomic micromotion orbit, the magnitude and direc- 
tion of the local magnetic field change in space and time, 
with the magnetic moment of the atom precessing around 
the direction of the field. In fact, given the Hamiltonian 
of Eq. a dynamical state \ty) satisfies the following 
center-of-mass equations of motion: 



d 2 K _ fj, 

dS fi 

dt s 



V[S ■ b(R)] - mgz 
(tf|sAB(x,t)|tf), 



(6) 



where R = Xx + Yy + Zz = |x|W) is the expectation 
value of the center of mass position, and S = S x x + S y y + 
S z z = (\I/|s|\I f ) is the expectation value of the atomic 
spin. 

Thus we are lead to a set of equations which are not 
closed. However if the quantum mechanical wave func- 
tion can be factorized, as to the lowest order approxi- 
mation, we may write (s A B(x, <)) = (s) A (B(x, t)) . 
Within this approximation the above system of equations 
becomes closed and assumes the form 



d 2 R 

'HP 

dS 
~dt 



V[S • b(R)] 
S A B(R, t) . 



mgz 



(7) 



The simplest periodic solutions generated by the Eqs. 
give a good estimate of the fast center-of-mass con- 
densate dynamics, i.e. the atomic micromotion. This 



motion was experimentally observed in a triaxial TOP 
trap in Ref. (Tj. This periodic solution is best viewed 
in a frame rotating with the bias magnetic field. In this 
frame the magnetic moment results to be aligned to the 
effective magnetic field and the atomic center of mass is 
at rest. Also, in order to have a stable motion, the spin 
must be tuned to the effective magnetic field, producing 
a confining potential energy. The stability along the z 
axis requires a zero force along this direction: this leads 
to 



S, 
s 



2b\n\ 



(8) 



which accounts for the fact that the component of the 
spin along z must be positive. The centrifugal force bal- 
ancing the gradient force yields the radius of the micro- 
motion 



\l*\b 



mujn 



(9) 



where £ and a are respectively the cosinus and sinus that 
the effective field B ms * forms with the z axis. The com- 
ponents of the effective field in the rotating frame, with 
horizontal component along the x axis, are 



b; 



br 



and Bl nst = -2bz + u T s/fj,. (10) 



In order to determine the z height of the periodic orbit 
under examination, the following parameters also useful 
in the subsequent analysis are required: 



luts 



B 



(11) 



p being the radius of the circle of death [l3| and <; is twice 
the amount the zero point of the quadrupole field shifts 
downwards for effect of the uniform fictitious magnetic 
field which appears in the rotating frame. In terms of 
these parameters the equilibrium height can be expressed 
by means of the following alignment relation 



-2b (z 



I) 



B + br 



which leads to 



2 sr^e 



B + br _ £ _ £ p 
1 2b ' 2 2 l a 



(12) 



(13) 



The above analysis suggest that the adiabatic approx- 
imation is more rigorous if we refer to B mst instead of 
the real magnetic field. Indeed, in the motion described 
above the adiabatic approximation is completely fulfilled 
with respect to this field. However, for more general so- 
lutions of the Eqs. Q) the behavior is different and we 
must also account for a spin component orthogonal to 
the local magnetic field. In any case, as it will be shown 
in the numerical simulations, the projection of the mag- 
netic moment along this field is a much better conserved 
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quantity with respect to that along the real field. In 
Refs. 0, HU the effect of the components perpendicular 
to the real magnetic field have been examined and it was 
found that a small misalignment with respect to this field 
gives rise to a Lorentz-type force. An additional electric- 
type force is originated by the time average of the fast 
oscillatory force induced by the spin precession. Both 
kind of forces, affecting the center-of mass dynamics, are 
geometric forces because they do not depend on the mag- 
nitude of the magnetic field but only on its orientation. 



III. FROM CLASSICAL TO QUANTUM 
ADIABATIC APPROXIMATION 

The gross features of the atomic confinement in a mag- 
netic trap are explained in terms of the adiabatic approx- 
imation. If this is fulfilled, an atom with electron mag- 
netic moment parallel to the local magnetic field experi- 
ences a confining potential given by U = |/zB|. Working 
in the rotating frame one replaces B with B ms *. Since 
in the following we will make use of an adiabatic-like so- 
lution in order to find non-adiabatic corrections, we will 
assume the spin to form an angle with B ms *, whence we 
will take the spin component in the direction of the in- 
stantaneous magnetic field to be sa = s • B inst /|B inat |. 
Thus, from Eq. we obtain the adiabatic Hamiltonian 
to be 



m 



p 

2m 



U a d , 



where 



U a d = mgz + <r|/x| (bx + Bq cosLOxt) 2 + 



(by + B sin L<j T t) 2 + (-2bz 



UJTS_^ 2 



1/2 



(15) 



For small displacements of the atoms from the equilib- 
rium position (0,0, h), the adiabatic potential U a d can 
be expanded in a power series of the displacement coor- 
dinates (x, y, £ = z — h), and up to the second order we 
have 



~m [cul r (x 2 + y 2 ) + ul z C 2 ] + U + U(t) . 



(16) 



Here the time-independent component is 

U Q = mgh + trp\/i\B , (17) 
and the time-dependent component is 

U(t,a) = ^^- [x cos u> T t + y sin u T t)- 

[(a; 2 - y 2 ) cos(2u T t) + 2xy su\(2uj T t)] . 



(18) 



The new adimensional constants here introduced are 

2h- cr 



V 



while the oscillation frequencies are 



u>o r — 



I a\p,\b(2 V 2 + 1) 



2mp(3 3 



Ua\n\b 



(19) 



Two time scales are involved in U. 



(2) 



ad ' 



y mpf3 3 

the slower one be- 
ing associated to the harmonic motion at the oscillation 
frequencies wo, r and u>o lZ , whereas the faster one is asso- 
ciated with the bias frequency u>t- The oscillating forces 
have vanishing time average over a period 2-k/lut- This 
can be substantiated by the following wave function fac- 
torization 



where 



*(x, t, a) = $(x, t, <r)£(x, t, a) 



£ = exp[— iw/h], 



(20) 



(21) 



and w — f Q dtU(t, a) describes the dominant effects of 
the oscillating potential. The time-scale separation al- 
lows us to consider $(x, t) as a slowly varying function 
of time j2|| . Notice the explicit dependence on the pa- 
rameter a. Substituting Eq. if^UI) into the Schrodinger 

(2) 

equation with the potential U ad , we get 



(14) ihdM*,t,a) = 



m[ojl r (x 2 + y 2 )+u 2 0z ( 2 



$(x,i,cr)- 



— Vw • V 

m 



2m 



ih 
2m 



(22) 



The above assumptions on the different timescales allow 
us to consider the coefficients of the oscillating terms at 
frequencies lot and 2lot as slowly varying ones. Indeed a 
time average over the short time 2tt / lot leads to 



ift<9 f $(x, t, a) = 



2m 



mg 




-m\LU 2 (x 2 + y 2 )+ui( 



where the irrelevant constant terms have been dropped. 
This equation displays a three-dimensional harmonic os- 
cillator structure whose frequences are 



1 



g (32t7 2 + 1) 



-,1/2 



l6£/3 3 pLo T {2tf 

1/2 



2/3 3 ptu T 



9 2 



(24) 
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in the xy-plane and along the z-direction, respectively. 
The equilibrium position along the z-axis is obtained 
by setting to zero the term multiplying £ in Eq. 1)23( 1 . 
Neglecting smaller contributions, the equilibrium atomic 
position is given by 



B. Local basis 



h{a) = s/2 - pt/(2^/a^e) . 



(25) 



In the limit a = 1 the classical solution of Eq. ((13)1 up to 
a small term containing w T 2 is matched. The expectation 
value of the particle momentum p on the state 1)20(1 is 



(p)= \l i \W° 3 _ ^(-sino^cosa^O). (26) 



LOT 



This result is only in part equivalent to the classical one 
of Eq. 1(13(1 because here the instantaneous position value 
(x) is always zero. 



IV. QUANTUM DYNAMICS 



It is useful to introduce a coordinate-dependent spin 
basis {|xm(x, t))} such that 



s-B R (x, t) 



|Xm(x, t)) — m|xm(x, t)) for 



j < m < j . 

The local basis vectors in which the z— axis coincides 
with the magnetic field in the same point can be given in 
terms of the angles $ and ip through the rotation operator 
M(ip, $) as follows 

\xMv))=M{<p,#)\j,k)=e-t* a 'e-i* a y\j,k). (31) 

With the total wave function expanded as |^ R ) = 
El=-j ^™( x , *) IXm(x, t)) , the Hamiltonian of Eq. ggj 
becomes 



if 



P , M nR/ x 
:- Vmgz B (x)s z + 

Zm s 

^-{2A.p + p(A)+A 2 }+V 



(32) 



A. Effective Hamiltonian 



with 



The quantum dynamics is more properly addressed 
by transforming the original equations into a spin ref- 
erence frame rotating at the bias field frequency. Thus 
the wave function \^ R ) in the rotating frame is writ- 
ten as |* R ) = R z (~LL> T t)\^) = exp(jrS z u T t)\^) , where 
R z ($) is the rotation around the z axis by an angle i9, 
and is the laboratory frame wave function. Then the 
Schrodinger equation for |>]/ R } becomes 



ihd t \ty R ) =ihd t R z {-u T t)\^>) = 
{H'(t)-huj TSz }\* R ) 



(27) 



where H'(t) 
and 



H'{t)-hw T s z = 



2m 



-mgz- -s-B R (x,i) , (28) 
s 



with B R (x, t) = [Bo + b{xcosu>Tt + ysmu)Tt)]x + 
b(ycosu>Tt — xsmuiTt)y + (—2bz + luts/h)z . B R is the 
magnetic field in the spin rotating frame, a constant bias 
field and a rotating inhomogeneous field. This effective 
magnetic field B R (x, t) = B K x + B R y + B K z identifies 
the position dependent angles and ip as 



i(B*y + (B*y 

■d = arctan — — r; , ip = arctan . (29) 



S R 



and B 



R _ 



{B*y + {Bfy + {Bfy. 



A = — (s z cos$ — s x sm-d)\7ip — Sj,W , 
V = — (s z cost) — s x smi3)d t p — s y dt$ ■ 



(33) 



Appendix ^ contains details useful to derive the func- 
tions in Eq. ((33(1 . In Eq. ((32|) A and V represent pseu- 
dopotentials connected to the Lorentz-like and electric- 
like kinds of forces introduced in Ref. [3(. The functions 
tp(x, t) and 7?(x, t) depend on the effective magnetic field 
geometry as stated in Eqs. (12911 . Berry's geometric terms 
appear in the quantum Hamiltonian through these an- 
gles. By direct computation the following relations are 
derived: 



R z (—u)Tt)H (t)R z (u>Tt) is time-dependent y<# = 



[B R A (B R A VB R )] . 



dt-d = 



(B R )V(S R ) 2 - (S R ) 2 
[B R A(B R A9 t B R )] z 

(B R ) 2 y / JW^) 2 - (_B R ) 2 



Vp = 



d t ip = 



[b r a vb r ] z 

(b r )2 - (Bjy 

[B R A9 t B R ] z 



(S R ) 5 



(34) 



These terms, invariant with the modulus of the magnetic 
field, are geometric fields depending only on the force 
lines of the magnetic field, i.e., the field geometry. Their 
explicit form is given in Appendix |B"1 

Notice that, starting from the Hamiltonian QJ, which 
is linear in the spin operators, the nonadiabatic terms 
give rise to a dynamics with a quadratic dependence in 
the spins. We may write i? R in a form that put this in 
more evidence 



H n = P 



2 

2m 



mgz 



^Vsi+^g^SjSk, (35) 
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where the indices (i,j,k) run on (x,y, z), and the spin 
coefficients are 



h x = — !— [ps'mdV ip + sindWipp] + sm-ddtf , 
2m 

h y = — — [pVtf + Vtfp] - d t d , 
2m 

h z — -—B R — [p cos -dVip + cost9V(/9pl — cos$<9 t <zj , 

s 2m 

g yy = J-|V^I 2 , 
2m 

^ „2 „q|Y7,„|2 



2m 



cos 2 z9|Vyf 



f = g zx 



2m 



sin^cos^lV^I 2 



(36) 



and g ]k — otherwise. Let us recall that all the operators 
hi,gjk are hcrmitian ones. 



V. BEYOND THE ADIABATIC 
APPROXIMATION 

A. Effective spin dynamics 

Since the exact solution of the Schrodinger Eq. I|27[) is 
an impracticable task, the spin dynamics will be taken 
into account in an effective way by resorting to a time- 
dependent variational principle (TDVP) [13. The TDVP 
procedure allows to reduce the system quantum dynamics 
to a semiclassical Hamiltonian form. This procedure was 
introduced for studying the low-lying collective states in 
nuclei |28| . but was later shown to provide a valid ap- 
proximation also for the one particle Schrodinger equa- 
tion. Within this procedure, whose details are shown in 
Appendix IO we choose a suitable trial state of the form 



V>(x,t) = exp[-iS(t)/h]V(x,t)\j,T(t)) 



(37) 



which will be subjected to the weaker form of the 
Schrodinger equation embodied into TDVP, i.e. Eq. (|C1|) 
in Appendix IO Here ^(x, t) and \ j, r(t)) take in account 
for the center-of-mass motion and spin dynamics, respec- 
tively. S(t) is an effective action for the spin variables. 
By carrying out the variational procedure on the trial 
wave function we derive the classical equations of motion 
for the expectation- values of the spin-operators Sj on the 
spin component of the dynamical trial state |j,r(f)). A 
key point in this variational procedure is the parametriza- 
tion of the spin variables in terms of coherent atomic 
states. These latter have the physical significance of an- 
gular momentum states produced by a classical source 
|28| . They depend on a complex parameter r and are 



defined as 



\j,r) 



3 

J m——j 



2j 

(j + m) 



1/2 



^+ m |j,m) 



where \j, n) are the spin basis with the quantization axis 
taken along the direction of the local field. These states, 
analogous to the coherent states of the electromagnetic 
field, are defined within a subspace determined by the 
angular momentum j. Within this sub-space each state, 
completely defined by the complex number r, is mapped 
onto the direction of a vector on a sphere by a projec- 
tive transformation j^- In our case this vector iden- 
tifies the orientation of a classical spin with respect to 
a local frame having the z axis along the local mag- 
netic field B R (x, t). This property can be understood 
by computing the expectation values of the spin compo- 
nents on \j, t) . By keeping in mind the parametrization 
t = -e-**' tan(i?'/2) we find 



= Ui T \ s x | j, t) = jh sin cos ip' , 
Sy = {j,T\s y \j,T) = j h sin tf' sin tp' , 
Sz = (j,T\s z \j,T) = jhcos'd' . 



(39) 



where and ip' are the angles between the classical 
spin and the local magnetic field. The detail of the 
spin dynamics derivation are contained in the Appendix 
IO Their ruling equations are generated by the classical 
Hamiltonian of Eq. (|C4|) . 

Let us focus on the center-of-mass motion described by 
the wave function ^/(x, t) as in Eq. 12U|) . The trapping 
potential obtained by the application of TDVP procedure 
as from Eq. I|C2|I can be expanded in a power series of 
the displacement coordinates around the trap center by 
keeping only the harmonic terms. As a matter of fact 
these terms depend on the quantity a(t), with a(t) given 
by 



s 



(40) 



that coincide with the definition introduced previously 
within the adiabatic approximation. Since the evolution 
of this quantity is much slower than the bias-frequency 
w t ( w z << w t)j its time dependence is maintained in 
the ruling equations even after averaging over the short 
time scale of the bias field time dependent terms as done 
in order to arrive to Eq. I|23|l . As a consequence the wave 
function solution of (|23|l can be written as 



$(x,i,cr) 



c„£ n (£)$ n (x, o) 



(41) 



where the vector index means n = (m, 11,2,113) along the 
three orthogonal directions, and the constants c n are de- 
termined by the atomic initial conditions. The functions 
$ n (x, a) are the eigenfunctions of the three-dimensional 
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harmonic oscillator with eigenvalues 



E n (a) =U 



1 / a/ib 
2m \Pujt 



+ 



(42) 



two r {a)(ni + n 2 + 1) + hjj z (a){n s + -) , 



and £ n (t) = exp[ij n (t) / h 
also a geometric phase 

7 n (£) =ih da I dx$, 

Ja(Q) 



i f Q dtE n ((j(t))/h] embodi 



(x, cr) 



ICS 



(43) 



The parameter <j(t) entering into the equations of motion 
for the atomic center of mass is actually a dynamical 
degree of freedom whose evolution is generated by the 
classical spin dynamics. Thus, the center of mass motion 
and the spin dynamics interact the one with the other and 
they must be simultaneously integrated. Let us stress 
that the dynamics we have just found, is the classical 
canonical counterpart of that one generated by the full 
quantum Hamiltonian written in Eq. (|35fl . 



B. Ground state configuration 

In order to find the non-adiabatic corrections to the 
ground state solution (J2DJ, we assume that this solution 
is well represented also if we keep the lowest order in 
the expression of the Hamiltonian of Eq. (|C4|) . i.e., the 
first order in (x,y,C) appearing in H R . This means that 
the Hamiltonian parameters of Eq. (|C4(I . the classical 
form of Eq. i|35|) , can be computed as an average on the 
adiabatic-like ground state solution of Eq. I|2(J|) of the 
approximated operators h z , g^ k . Therefore at the lowest 
order of approximation the terms in Eq. I|3bj) result 

h x ~ -[— sm(ojTt)p x + cos(cuTt)py] , 



h v 



mp/3 
1 



mp(3' 



-{r)[cos(u) T t)p x + sm(w T t)p v } + 2p z } , 



V 



s mpp 
1 



sm(u> T t)p x + cos(u>Tt)Py\ , 



ff 



2mp 2 f3 2 ' 
3 + P 2 
2mp 2 P 4 ' 

2mp 2 (3 2 
V 



2mp 2 (3 2 ' 

Recalling the expectation value of the momentum p given 
by Eq. up to the order 1/p we have 



2(3 u>i r 



(l+V 2 ) U T ' 

pB _ 2r]/3 ujl r 
(l + V 2 ) ut 



and otherwise. The corresponding classical spin Hamil- 
tonian is, apart a constant, 

whose equations of motion results 
S x = -\hfo))Sy , 



S y — —{h*ty)S z + (h z ^)S x 
S z = {h* -))Sy ■ 



(44) 



By setting Si — with (i = x,y,z) we determine the 
ground state configuration 



(SlS° y , <S°) = (± 



hj 



s 

V 1 + ( /i (0)) 2 /(^ )) 2 ' 



= ,0, 



(45) 



The S° solution, with the spin aligned to the local mag- 
netic field, provides a correction to the adiabatic approx- 
imation discussed before. 



C. Effective classical dynamics 

The center of mass motion is described by the wave- 
function ^(x,t) = 3>(x,t, 0")£(x, i, it) introduced in lET7)l . 
The exponential factor £(x, t, a) has been defined in 
Eq. (|21(l and $(x, t, a) satisfies the time-dependent 
Schrodinger equation l|23|) for the 3-dimensional har- 
monic oscillator which Hamiltonian is H = p 2 /(2m) + 
C//j(x, a(i)). In terms of the frequencies (|24|l and 
of the equilibrium atomic position l)25[l. the time- 
dependent harmonic potential has the form t//j(x, a(t)) = 
m[w 2 {a{t))(x 2 + y 2 ) + iJ 2 z (a(t))(z ~ h(a(t))) 2 }/2. 

Upon introducing the center of mass position R = 
|x| SEf) and momentum P = (^Ipl^), the following clas- 
sical equations of motion are easily derived 



dR 

~dt 
dV 

~dt 



-V R U h (R,a(t)) 



dAP 

dt 



(46) 



where AP = J dx\$(x, i)| 2 Vw(x, t), and a(t) is defined 
by Eq. (gOJ). Now we introduce a further factorization 
of the kind (*|C(x,t)p|*) w 0(R,i)P, where 0(x,£)p 
stands for the first three among the operators appearing 
into Eqs. I|36(l expanded in a power series of x and p up 
to the second order |29j. Then by considering only the 
linear terms in spin variables appearing into the classical 
spin Hamiltonian l|C4f) we can write 



(47) 
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where the time-dependent coefficients are implicit in the 
center of mass wave function ^(x, i) as given in JCSJ. 
Thus in Appendix |0 we derive the following equation of 
motion for the classical spin 



is well conserved while the spin projection along the real 
magnetic field B evidences time dependent oscillations, 
as already stated by Ref. . 



dS 
~dt 



= B(t)AS. 



(48) 



where the magnetic field B is the sum of the real one plus 
some fictitious terms having originated from the geomet- 
ric forces with components 



B x (t) 



{B*A[(£.V + 9 t )B*]} ; 



By(t) = - 



{B R A[B* A [(£.V + a t )B*]]} z 



(B*)\(Bfy + {Bfy 



E(f)= M pR {B R A [(£ ■ V + <9 t ) B R ] } z 

Z{) s B R [(B, R ) a + (B») s ] 

(49) 

Thus, the two equations systems (|46f) and (|48|l form a 
closed system to be simultaneously integrated. 



VI. NUMERICAL SIMULATIONS 

We have numerically integrated the set of equations 
(|46|1 and l)48|) by means of a Runge-Kutta algorithm. 
The set of parameters chosen, i.e. Bq = 4. 10~ 4 T, 
b = 0.18T/m, and lot — 2n 10 4 s _1 , correspond to those 
used in TOP experiments exploring the rubidium micro- 
motion 0, 0| . The simulations allowed to recover the 
atomic micromotion, representing periodic closed orbits. 
The micromotion was investigated through the classical 
equations of motion J7J and also through the improved 
system of equations of Eqs. I|4tj|) and (|48|l . Similar results 
were obtained for the center of mass motion. In both ap- 
proaches we observed a strong dependence on the initial 
conditions, that for the classical center of mass variables 
are given by Eqs. i|25|l and (|26|l . For the spin variable 
the classical condition corresponds to the spin aligned 
along the local B ms * magnetic field, while the quantum 
mechanical solution requires the atom to be in an eigen- 
state of the spin operator along the local magnetic field. 
A modification of the initial conditions from those re- 
quired for the atomic micromotion, for instance a shift 
of 100 /im along the 2— axis, produced the open trajec- 
tories shown in Fig. [21 We noticed also a strong depen- 
dence on the initial condition for the atomic spin. We 
also verified that within the parameters used here which 
approximately match those corresponding to the experi- 
mental set up of Ref. [llj] , the correction to the adiabatic 
approximation expressed by Eq. 145fl are not quite rel- 
evant. We verified numerically that the spin projection 
along the effective magnetic field B ms * given by Eq. (|10|l , 




5303-0.003 



FIG. 2: Plot of a 3-D trajectory originated by initial condi- 
tions close to those corresponding to the micromotion. The 
numbers on the axes express the coordinates in fj,m. This 
trajectory corresponds to 5 times the 2-k/ujt period. It is 
not closed, while a stable micromotion orbit corresponds to a 
closed motion. The number on the axes express the coordi- 
nates in fim. 

We have explored a different region of parameter val- 
ues, where we expect the adiabatic approximation to 
break down. While the adiabaticity is certainly not ful- 
filled if lu z , to r ~ wj, , a further source of failure for this ap- 
proximation rests in the intensity of the geometric mag- 
netic fields of Eqs. (|49|l being of the same order of magni- 
tude of the applied real ones. This occurs, for example, 
at bias field intensities of the order of Bq = 2. 10~ 7 T. 
Notice that for this weak rotating RF field, the trap oscil- 
lation frequencies of Eq. (|19|l are large enough to sustain 
the atoms against gravity. However the radius p of the 
circle of death of Eq. (|11J) becomes comparable to the 
radius r of the micromotion orbit. Fig. 03 (a), (b) and (c) 
show the components of B as a function of time obtained 
with initial conditions very close to those of a micromo- 
tion orbit. For this set of parameters it is interesting to 
make a comparison between the atom dynamics gener- 
ated by the classical equation of motion and the effective 
improved equations of motion. The latter give rise to a 
stable motion, shown in Fig. |3fd) traced by integrating 
the effective equations l|46l) and (|48J) . Under the same 
initial conditions the classical equations of motion gen- 
erate an unstable trajectory with the atoms conserving 
initially a constant height z = 0.73/xm, and then after 
several milliseconds escaping from the trap. The more 
stable character of the effective equations solution in re- 
spect to the classical ones is made evident by comparing 
the spin projection along B ms * in both cases. By numer- 
ically integrating the effective equations (|46|l and 14811 . 
we found oscillations of S ■ B ms */I<S>||B™ S *| near the sta- 
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ble value 1. On the contrary, by integrating the classical 
equations Q, we found spin flip that causes the conden- 
sate escape from the trap. 




time (ms) 



FIG. 3: In (a), (b) and (c) plots of the B x , B y , and B z geo- 
metric fields, in units of 1. x 10~ 9 T, as a function of time in 
ms for atomic motion within a TOP trap with a RF bias field 
of 200 nT. In(d) the z— position, in /jm, of the atomic center 
of mass traced by integrating the effective equations l|46|l . In- 
stead the integration of the classical equations of motion 
displays an unstable trajectory. 

The important role played by the terms originated by 
the non adiabatic approximation appears very clearly 
when we compared the atomic equilibrium within the 
TOP z cq as derived by the classical solution to that 
predicted by the effective Eqs. and |gSJ|. That 

comparison is shown in Fig. 0] for a fixed quadrupole 
field b = 0.18T/m and a RF rotating field Bo between 
1.x 1(T 4 T and 2. x 1(T 7 T. At large values of B the values 
z eq predicted by the classical solution and the improved 
one coincide. Instead the two values are different at small 
values of B because the two solutions predict different 
equilibrium positions. Finally, for Bo < 5. x 10~ 7 T, clas- 
sically the atoms are not suspended against gravity, while 
the effective equations predict a stable equilibrium posi- 
tion. 



VII. CONCLUSIONS 

In this paper we have analyzed the motion of neutral 
atoms within TOP magnetic traps. We have consid- 
ered the approximate classical equations of motion de- 
scribing such system and revisited the fast degrees of 
freedom motion at the forcing frequency lot known as 




FIG. 4: Comparison between the equilibrium positions z eq 
of rubidium atoms within the TOP trap versus the B D RF 
field as predicted by the classical solution and by the effective 
equations 146H and 14811 at a given value of the quadrupole 
field gradient b = O.I8T/m. For B < 5. x I0~ 7 T, the classical 
model does not leads to a stable orbit. 



atomic micromotion. In such a motion non-adiabatic ef- 
fects and geometric fields are absent. Within the scenario 
of adiabatic approximation, we have analyzed the atomic 
quantum motion by taking into account the lowest fre- 
quencies terms embodied into the time-dependent poten- 
tial. Within this approximation, the center-of-mass mo- 
tion resulting harmonic, we have calculated the harmonic 
trapping frequences and recovered the quantum counter- 
part of the classical micromotion. Addressing the prob- 
lem within a pure quantum context, as a consequence of 
the presence of inhomogeneous magnetic fields geometric 
magnetic fields appear quite naturally. These geometric 
fields are responsible for a misalignment of the atomic 
spin with respect to the local magnetic field and then for 
non-adiabatic effects. Within this framework, we have 
derived an effective classical dynamics in which these geo- 
metric fields are explicitly embodied. The atomic motion 
results by the coupling of a quantum harmonic motion, 
governing the center of mass, and an effective nonlinear 
spin dynamics driven by both the local magnetic field 
and the geometric ones. 

The numerical simulations performed for the param- 
eters of standard experimental set-ups have shown that 
the adiabatic approximation is well suited. On the other 
hand, by reducing the intensity of the bias field, non- 
adiabatic effects show up, because the geometric field 
becomes more relevant and cause misalignment of the 
spins around B ms *. 

Another relevant facet concerns the sensitivity to the 
initial condition of the trap equations. The projection 
of the atomic magnetic moment on the local field B mst 
is a conserved quantity, an adiabatic integral. The ini- 
tial condition of S • B mst chosen for a given simulation 
identifies a given dynamical evolution, its value remain- 
ing conserved during atomic motion. On the other hand, 
the equilibrium height is determined by the value of this 
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quantity [see Eq. flUM, so that the projection of the 
atomic magnetic moment on the local field determines 
the cloud equilibrium height. Now, it is a fact that, for 
a given geometry of magnetic fields within a trap, such 
a height results to be independent on the initial prepa- 
ration condition. Thus, since the dynamical equations 
do not select by itself any special value of this quantity, 
one could argue, in order to explain some experimental 
features of the BEC clouds 0, , that a possible mod- 
ification of the actual spin projection on B ms< might by 
involved in some steps of the preparation of the relative 
Bose-Einstein condensate. Further attention should also 
be devoted to the non-linear interaction within the con- 
densate. 



VIII. ACKNOWLEDGMENTS 



This work has progressed through the constant collab- 
oration of several members of the BEC group in Pisa, M. 
Anderlini, D. Ciampini. R. Mannella, O. Morsch, J.H. 
Miiller, and through useful discussions with M. Mintchev. 



APPENDIX B: GEOMETRIC TERMS 

For a cylindric TOP trap with magnetic field as defined 
in (0), the geometric fields result 



Vi? 
Vp 



(5 R )VW) 2 + (B?) 2 
b 2 



Bf b(x + pcosuirt) 
Bf b(y + psmu) T t) 
2[(B R ) 2 + (Bf) 2 } 



(BR)* + {BfY 
u; T b 2 Bf 



-y — psin tuxt 
x + p cos LOrt 




d t (p = LO T 

Ad = 



(B*) 2 J(B*) 2 + (Bf) 2 



: p (y cos tort — x sin urt) . 



(5«) 2 + (B^ 



- 1 



Vy>W = , Aip = . 



1 + 6 



(B R ) a 



(Bl) 



APPENDIX C: TDVP APPROACH 



APPENDIX A: ABOUT THE SU(2) ALGEBRA 



The su(2) algebra is defined starting from the angular 
momentum generators s x ,s y ,s z and their commutation 
relations. Using the standard definitions for the raising 
and lowering operators, we derive the following relations 
for the derivatives of the operator M(<f>, 9) defined in Eq. 



M\p,d)d ip M{p,d) =-t:(s* cost? - S:E sin tf), 

M^(ip, &)d#M(<p, i?) = — ^s y , 
M \p, d)d 2 w M{v, 0) = - ^(s z cosi? - s x sintf) 2 , 
M\p,d)d 2 § M{p,$)=-^s 2 y , 

M f (tp, 0)dfoM (<p, i?) = - (s z cos tf-s x sin $)s y . 

M'(if,i})s x M(ip,i}) —SxCoadcostp — s y sinip+ 

s z sin & cos tp , 
M'(tp, <&)syM((p, $) =s x cosz9sint/3 + s y cos<£>+ 

s z sin i? sin tp . 



In our contest the TDVP method structures the dy- 
namical quantum-state describing the atomic motion, 
in terms of a trial state written, as in Eq. I|37(l . as 
the product of a time-dependent phase factor e ~ lS ( t )/ h 
times a spatial- and time-dependent scalar wave func- 
tion ^(x, t) times a time-dependent spinor \j,r(t)). The 
time-dependent trial state ^(x, t) is to be found in a 
self-consisting way. By imposing the weaker form of the 
Schrodinger equation 

J dx (x, t) [ihd t - H R ] ip(x, t) = , (CI) 

we get the effective action S(t) 

S(t) = f dt 
Jo 

J dx(7,r(t)|*t(x,t) [ihd t -H*] *(x,t)tf,T(t)>. 

By splitting the Hamiltonian as H R = H° d +AH a , where 
H° d has been introduced in JT3J, we get 



S(t) = J dt J dx*t(x,t) [ihdt - H° d ] *(x,t)- 



dt I dx(j,T(t)\ [ih\^(x,t)\ 2 d t - 

*t( x , t )AiTtt(x,t)] \j, T (t)}. 



(C2) 



(Al) 



By expanding the Hamiltonian H£ d in a power series of 
the displacement coordinates (x,y,£ = z — h), we get 
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the harmonic Hamiltonian with the time-dependent po- 
tential of Eq. (|16|) . Thus, by introducing the structure 
^(x, t) of Eq. Il2()[l for the atomic wave function we find 
the Schrodinger Eq. (|22[) . After taking the time aver- 
age of the latter equation on a short time 2it/u;t, we get 
the harmonic problem (|23|) of which the general solution 
$(x, t) is known. Thus, the first term in the r.h.s. in Eq. 
J2I) vanishes, and we obtain 



After boring algebra we get the equations of motion 



(It 
~dl 



(1 



|2\2 



2ihj 



drH(T,f), 



S(t)= / dt(j,T(t)\ im~(AH( 2) ) \j,r(t)), 
Jo L J 

where (Aff? 2) ) = / dx [*t( X) t)AH^(x, t)} and Aff^ 
is obtained by expanding AH a in a power series in x and 
p up to the second order. Explicitly we have 



where 7i(r, r) = (r| (AH^) |r) and r is the complex con- 
jugate of t. These tricky equations can be cast in more 
familiar form by writing them in terms of the classical 
spin components S x ,S y ,S z introduced in the text. The 
dynamics of these spin components, apart terms irrel- 
evant for the equations of motion, is generated by the 
following Hamiltonian: 



i jk 



Sk 



(2) 



with the obvious meaning of index (2). 

The Hamiltonian [AHZ^) is built of su(2) algebra gen- 
erators acting on the time-dependent spin vector \ j, r(t)}. 
We choose for | j, r(i)) the components of the su(2) atomic 
coherent state \r(t)) of j representation. As a support 
for this choice let us remember that, if the Hamilto- 
nian (AH^) was a closed dynamical algebra, i.e. a lin- 
ear combination of the su(2) generators, the solutions of 
the Schrodinger equation given by these coherent states 
would be exact J27J • The equation of motion for the la- 
bel r(t) (and its complex conjugate), that involves the 
dynamical evolution of the state ip(x,t), is obtained by 
stationarizing the effective action S(t) 



(C4) 



H(S x ,S y ,S z ) = y^ J (h} 2 ))S i + 

i 

where the time-dependent coefficients are derived by the 
center of mass wave function \I/(x, t) as 



(h\ 2) ) = y'dx[* t (x,i)/i| 2) *(x,t)] foii = x,y, 



(5(2)) = / dx [* T ( x ^)ff( 2 )*( x ^)] f or i,j = x,y,z. 



SS 




ih(j,r(t)\ 



dr dr 



\j,r(t))- 



(j,T(t)\(AH^)\j,r(t)) 



= 0. 



(C5) 



the operators h l and g % i have been given in l|36(l . There- 
fore we obtain a canonical system for the classical spin 
components whose dynamics is generated by the effec- 
tive Hamiltonian Ti., endowed with the Poisson brackets 
{<S x ,<S y } = S z and cyclic permutations. 
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